Simulating the coupling of angular momenta in distant matter qubits 
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We present a mathematical proof of the algorithm allowing to generate all - symmetric and 
non-symmetric - total angular momentum eigenstates in remote matter qubits by projective mea- 
surements, proposed in Maser et al. [Phys. Rev. A 79, 033833 (2009)]. By deriving a recursion 
formula for the algorithm we show that the generated states are equal to the total angular mo- 
mentum eigenstates obtained via the usual quantum mechanical coupling of angular momenta. In 
this way we demonstrate that the algorithm is able to simulate the coupling of N spin-1/2 systems, 
and to implement the required Clebsch-Gordan coefficients, even though the particles never directly 
interact with each other. 

PACS numbers: 42.50.Dv,42.50.Tx,37.10.-x,03.67.-a 
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I. INTRODUCTION 

The generation of multipartite entangled states relies 
typically on a direct interaction between the particles, in- 
duced via nonlinear effects [5J [5] , collisions [2j , Coulomb 
coupling [SHI] or atom-photon interfaces [8j. However, 
recently it was shown that it is also possible to entan- 
gle distant particles that never directly couple to each 
other. Examples are entanglement swapping of entan- 
gled photons [9] produced via Spontaneous Parametric 
Down Conversion (SPDC) [2] [3] or entanglement of re- 
mote atoms generated via projective measurements of 
emitted photons [TOHIB"] . By use of the latter technique, 
distant clouds of atoms have been entangled so far [TOT - 
IT2"] . but it was also possible to demonstrate the entan- 
glement of distant single atoms, e. g., two ions stored in 
different ion traps one meter apart 13]. 

Recently, an algorithm was proposed, based on the 
technique of projective measurements, which allows the 
generation of all - symmetric and non-symmetric - to- 
tal angular momentum eigenstates in N distant matter 
qubits [16] , most of them being entangled states [17l [18] . 
The algorithm was motivated by comparison with the 
familiar construction in quantum mechanics of the eigen- 
states of a total angular momentum, obtained via cou- 
pling of N spin-1/2 systems, and tested numerically for 
up to N — 8 qubits. 

In this paper, we present a rigorous mathematical 
proof of this algorithm by the technique of induction. 
For the proof, we derive two recursion formulas, one for 
the algorithm presented in |16j and the other one for the 
familiar quantum mechanical procedure to generate the 
eigenstates of a total angular momentum of an A-qubit 
compound. By complete induction it is then shown that 



'Electronic address: joachim.vonzanthier@physik.uni-erlangen.de 
URL: http : / /www. optik .uni- erlangen. de/j vz/ 



both recursion formulas lead to the same result. In this 
way, we prove that the algorithm presented in [16] al- 
lows to generate all symmetric and non-symmetric total 
angular momentum eigenstates, and to implement the 
required Clebsch-Gordan coefficients, for any number A 
of spin-1/2 particles even though the particles never di- 
rectly interact with each other. The algorithm is thus 
able to simulate the coupling of angular momenta even 
of non-interacting qubits. 

The paper is organized as follows: in Sees. II and III 
we recall the algorithm presented in |16| and formulate 
it in a rigorous mathematical way. In Sec. IV we exem- 
plify the algorithm for a specific case and, in Sec. V, we 
rewrite it in form of a recursion relation. In Sect. VI we 
recapitulate the well-known quantum mechanical rules 
to generate the eigenstates of a total angular momen- 
tum from the eigenstates of the angular momenta to be 
coupled and rewrite it for the special case of A spin-1/2 
particles also in form of a recursion formula. In Sec. VII 
the two recursion formulas of Sec. V and VI are then 
compared and proven to be identical by the technique of 
induction. In Sec. VIII we finally conclude. 



II. ALGORITHM FOR THE GENERATION OF 
TOTAL ANGULAR MOMENTUM EIGENSTATES 
OF N SPIN-1/2 PARTICLES 

For A spin-1/2 particles, the total angular momentum 
eigenstates, defined as simultaneous eigenstates of the 
square of the total spin operator and its z-component 
Sn,zi are commonly denoted by \Stf', m^), with the cor- 
responding eigenvalues Sn(Sn + l)^ 2 and m^ti, respec- 
tively |17[I18|. However, the denomination \Sn', tun) may 
in general characterize more than one quantum state. We 
will thus specify more precisely an A qubit total angular 
momentum eigenstate by incorporating in its notation 
its coupling history, i.e., by writing the quantum num- 
bers Si, S2, Sif-i of the total spins Si, S2, Sjv in 
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addition to Sn and mjv- According to the quantum me- 
chanical rules of adding angular momenta, a single qubit 
state has Si = |, a two qubit system can either have 
S2 = or S2 = 1, a three qubit system S3 = \ or S3 = §, 
and so on. Including the coupling history we thus get the 
notation \Si, S%, Sn; m N ) for an eigenstate of the total 
spin S jv which describes a particular angular momentum 
eigenstate unambiguously [IB] . 

To generate all total angular momentum eigenstates in 
non-interacting matter qubits, Maser et al. considered 
in 16 A uncorrclatcd single photon emitters with a A- 
configuration, e.g., trapped neutral atoms, trapped ions, 
quantum dots, or any other equivalent physical system 
with access to a similar behavior. In the beginning each 
emitter is assumed to be in the excited state |e). Sub- 
sequently, due to spontaneous decay, the emitters leave 
the excited state |e) and populate the two stable ground 
states |— ) and |+) under the emission of a <x+ and er~ 
polarized photon, respectively; hereby, the two states |±) 
encode the qubit, with eigenvalues ±h/2 for the operator 
Si jZ . The A photons are either registered in the far field 
by A single photon detectors or collected and transmit- 
ted via optical single mode fibers towards the A detectors 
where it is assumed that each detector registers exactly 
one photon. Hereby, each photon can possibly travel via 
A paths towards the A detectors so that after a success- 
ful detection event it is impossible to determine along 
which way a recorded photon has travelled. This causes 
multi photon intereference effects of Ath order which can 
be exploited to generate a large variety of entangled multi 
qubit states [T3HIB] . 

Using this scheme and assuming that optical fibers for 
the transmission of the photons towards the A detectors 
are used, each accumulating an optical phase 8j k = 27m 
(with neZ) between emitter k and detector j, the al- 
gorithm for the generation of any of the 2^ symmetric 
and nonsymmetric total angular momentum eigenstates 
I Si, S2, S3, Sn', mjsr) spanning the Hilbert space of the 
A-qubit compound can be formulated in the following 
way [15] : 

1. set up y + rriN — itin) detectors with cr~ (<x+) 
polarized filters in front. Hereby, we connect the 
first emitter with optical fibers to all A detectors. 

2. check for each emitter k beginning with k = 2 
whether Sk > Sk-i or Sk < Sk-i- If 

a. Sk > Sk-i, we have to connect the emitter 
with optical fibers to all detectors except those 
which are mentioned in case [b.] below. 

b. Sk < Sk-i, we have to connect the emitter 
with optical fibers to one detector with a <t~ 



polarizer and to one with a <t + polarizer where 
the optical fiber leading to the <x~ polarizer 
should induce an additional relative optical 
phase shift of tt. Those two detectors should 
not be linked to any other subsequent emitter. 
If one wants to create a particular total angular mo- 
mentum eigenstate \Si, S2, S3, Sn', ttin), the setup to 
generate this state is unambiguously determined by the 
quantum numbers S\, S2, S3, Sn obtained by succes- 
sively coupling the spins of A spin-1/2 particles, and 
the magnetic quantum number m^v of the total spin Sn- 
Hereby, according to the rules of the algorithm presented 
above, mjv determines the number of cr~ and cr + polar- 
ising filters used in the setup. 



III. MATHEMATICAL FORMULATION OF 
THE ALGORITHM 

In order to prove the algorithm presented in [16] and 
recalled in Sect. II, we have to formulate it in a more 
mathematical way. To this end, we note that if emitter 
k emits a photon and we register this photon by any of 
the A detectors j (j,k € {1,..., A}), which is however 
not precisely known, equipped with a polarising filter in 
front with orientation ej = ctj<y + + fij<J~ , the state of 
our atomic system is projected by applying the following 
projection operator Pf. [TH ITS] : 
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( aj \~) k {e\+ P 3 \+) k (e\)®\\) 1 {0\. (1) 



Hereby, the phase factor e lSj - k accounts for the optical 
phase dj^k accumulated on the way from emitter k to 
detector j and the operator |±) fc (e\ projects the fcth 
atom from the excited state |e) to the ground state |±), 
while leaving the state of the other atoms unchanged. In 
Eq. (T| the operator 1 1)^- (0| ensures that every detector 
registers exactly one photon. While this is not relevant 
when considering only the photon emitted by a single 
atom it becomes important when each of several atoms 
scatters one photon. Since it is assumed that initially 
all atoms are in the excited state and all of the photonic 
modes are unpopulated, the initial state \ipi,N) of the 
atoms and the modes is given by 



I0...0), 



(2) 



After recording all A photons at the A detectors, the 
final state 1 4>f) of the system can be written in the fol- 
lowing form 
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N 1 

H kUN ~k)\ ^ ^ Pa(k)aa(k+l)---a a(N) Xa(l)r(l)---Xa(N)r(N) 
k=0 '' rGSjv ctGS.y 

(|+) T(1) (e| ® • •• ® |+) T(fe) (e| ® |-> T(fc+1) (e| ® ■■ ■ <g> |-) t(jv) (e|) l^.jv) 

V 

projection operator part 



with e i5 -o).-( fc ) = X^oxfe) and j,fc e {1,...,JV}. Hereby, 
the second and third sum run over all permutations er, r 
of the set {1, ...,N}, which can be mathematically de- 
noted by the elements of the symmetric group £/y [IS]- 
The summation index fc of the first sum controls the 
number of /3— and a— coefficients, i.e., the number of 
projection operators |+) (e| and |— ) (e|. Since we sum 
over all permutations of the projection operator for ev- 
ery k G {0,...,N}, we have to compensate the multi- 
ple counts in the sum J2 T e^ N the factor k \^ N _ k y ■ 
there are fc! combinations to choose the same set of pro- 
jection operators |+) (e| and (N — fc)! possibilities to 
choose the same set of projection operators |— ) (e|. The 
sum EctgSjv P<r{l) ■ ' ' Ar(fc)a<r(/c+i) • • ■ u a(N) must not be 
compensated for multiple counts, because these multiple 
counts are already included in the product Ph ■•■Pi, as 
each Pk is a sum over all a's and /3's (cf. Eq. |IJ). 

In general, the phase factors Xjk may take any complex 
value with norm between zero and one in order to allow 
for an attenuation of the light fields propagating from 
emitter fc to detector j. However, in the presented algo- 
rithm Xjk needs only to take the values 1 ("normal", i.e., 
not attenuated connection between emitter and detector 
with Sj y k = 27m, where n £ Z), (no connection between 
emitter and detector, meaning an attenuation of 100%) 
and —1 (not attenuated connection between emitter and 
detector with an additional phase shift of 7r with respect 
to a normal connection). 



I 

Note that we omitted in Eq. (|3]) the operators (0|. 
They were necessesary to derive Eq. ^ but would in 
the final form of Eq. (|3| only add a global factor of 

IljLo l-Oj which d° es n °t influence the final atomic 
state. Therefore, we chose to omit this factor for the sake 
of simplicity and accordingly change the initial state to 

foM) = |e...e). (4) 



IV. EXEMPLIFICATION OF THE ALGORITHM 

As an example for the operation of the algorithm, 
let us generate the three qubit state |; 5), where 

we have Si = 1/2, S2 = 1, S3 — 1/2 and 7713 = 
1/2. This state can be written in the decoupled basis 
as ^(2 |+ + -) - |+ - +) - |- + +)) [17 . Following 
the prescriptions of the algorithm of Sect. II we need 
Y + rriN — § + g = 2 cr — polarising filters in front of 
any two out of three detectors. As the index k of the first 
sum in Eq. ^ corresponds to the number of polarisers 
with er - — orientation, the only non-zero contribution to 
the sum comes from the summand with k = 2. If we use 
the well-known cycle notation from group theory |19j for 
the elements of the symmetric group S3, we can expand 
Eq. ([3| in the following way: 
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/ 

/3l/3 2 a-3XllX22X33 + ^2/?lQ!3X2lXl2X33 + ^3^2^1X31X22X13 

S v ' S - ' V - ' 

y a—id.r—id a — (12),r—id a — (13), r— id 

+ ^l^ 3 0;2XllX32X23 + /?2/33aiX2lX32Xl3 + ^3/3ia2X3lXl2X23 |+ + ~) 



<r=(23), T=id 



ff=(123),T=«2 



cr=(132),T=id 



+ - 



1 



/ 

^1^2Q!3Xl2X2lX33 + ^2/3l<33X22XllX33 + ^3^2«lX32X2lXl3 

v v ' v s, ' " v ' 

\ cr=id,T=(12) ct=(12),t=(12) cr=(13),r=(12) 

+ ^l^ 3 a 2 Xl2X3lX23 + /32^3aiX22X3lXl3 + /53/5lQ!2X32XllX23 |+ + ~) 



<t=(23),t=(12) 



cr=(123),T=(12) 



tr=(132),r=(12) 



+ - 



1 



/ 



/3l/32«3Xl3X22X31 + /32/5l«3X23Xl2X31 + /33/32aiX33X22X 
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\ a=id,T={13) ct=(12),t=(13) <t=(13),t=(13) 

+ /3l/?3a2Xl3X32X21 + /32/33«lX23X32Xll + /33/3l02X33Xl2X21 



(J=(23),t=(13) 



= (123),t=(13) 



<t=(132),t=(13) 



/ 

/3l/320-3XllX23X32 + /52/5l«3X2lXl3X32 + /33/32aiX3lX23Xl2 
v v ' v v ' v v ' 

\ a=id,T=(23) <t=(12),t=(23) ct=(13),t=(23) 



(5) 



/31/3 3 a-2xllx33x22 - 

S v ' 


|-/32/3 3 OlX2lX33Xl2 - 


f ^3/3l02X3lXl3X22 




<t=(23),t=(23) 


ct=(123),t=(23) 


' 'J 

<t=(132),t=(23) / 




/3l/?2Q:3Xl2X23X31 4 

v v ' 


-^2/3l03X22Xl3X31 4 


-/3 3 /3 2 aiX32X23Xll 


CT =id,T=(123) 


ct=(12) : t=(123) 


ct=(13),t=(123) 




^1/3302X12X33X21 - 

" v ' 


l-/3 2 /3 3 a;iX22X33Xii - 


f /33/3ia2X32Xl3X21 




ct=(23),t=(123) 


cr=(123),T=(123) 


ct=(132),t=(123) / 




/3l/32«3Xl3X2lX32 4 
s * ' 


-/32/3l03X23XllX32 ~\ 


-/3 3 /3 2 aiX33X2iXi2 

S v ' 


CT =id,T=(132) 


ct=(12) : t=(132) 


ct=(13),t=(132) 










/3l/3 3 02Xl3X3lX22 " 

V v ' 


)"/32/33aiX23X3lXl2 - 


f /3 3 /3iQ!2X33XllX22 




ct=(23),t=(132) 


cr=(123),T=(132) 


£7=(132),t=(132) / 





According to the recipe of the algorithm we could set up 
the polarization filters and connect the emitters to the 
detectors for example according to the scheme depicted 



in Fig.[TJ 

From the figure, we can now determine the required 
values for the coefficients oij , j3j and the phase factors 
Xjk appearing in Eq. (15]). There are two cr~— polarising 
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the state we want to create we have to include these 
quantum numbers to the state on the left hand side of 
Eq. (|3j). This allows us to know exactly the particular 
values we have to choose for otj, j3j and Xjk in order 
to implement the targeted state |Si, Sn] itin) by use 
of the algorithm. Later on, we will also use a notation 



of the form P, 



N ' 



Pi \4>i, 



N, 



, which means that the 

sum in Eq. (I3J) over the index k reduces to the summand 
with k = x. 



FIG. 1: Setup for the generation of the state ~, 1, ~; 



RECURSIVITY OF THE ALGORITHM 



filters, one in front of detector 1 and the other one in 
front of detector 3, thus we have 



a 3 =0 A ft = & 



1. 



(6) 



In front of detector 2, a <x+ — polarising filter is positioned 
and therefore we have to set 



a 2 



1 A h 



0. 



(7) 



A normal connection between an emitter k and a detector 
j has an optical phase factor Xjk — 1 (depicted in Figfl] 
by a black line). A missing link indicates an attenuation 
of 100% and thus Xjk — 0. Last, the optical fiber loop 
(labeled in Figjljwith n) indicates an additional optical 
phase shift of 7r and hence Xjk = — 1- In short we thus 
get: 



XU = X21 =X31 = Xl2 = X22 = X32 = X23 = 1 

X33 = -1 A Xl3 = 



(8) 



Hence, we can calculate the final state by inserting 
Eq. I^-Q into Eq. what leads to: 



P3P2P1 1^,3) 



|i 1 

1 2 2 ' 2 / 



(9) 

corresponding - up to a normalization constant - to the 
desired state 1 1 , 1 , | ; § } • 

Note that we introduced in Eq. Q a new notation 

when writing the state as Pn ■ ■ ■ P\ \iP%,n) 

\Si,...,S N ;m N ) 

Hereby, we took into account that the quantum numbers 
S%, S2, ■-, Sn, labeling the state |<Si, Sn', win), deter- 
mine unambiguously the values for otj, /3j and Xjk (with 
j,k E {!,... ,N}). Thus, in order to describe exactly 



Maser et al. have shown in [16] that the generation of 
the state 1 1, 1, |; |) following the recipe of the algorithm 
is in fact obtained in two steps. First, the states 1 1, 1; l) 
and ||,1;0) are generated and then, in a second step, 
these two states are superposed by use of another qubit 
|| ; ±|), just like quantum mechanically coupling a spin- 
1 system with another spin-1/2 system |16j . 

Inspired by this procedure, we will try to retrieve a gen- 
eral recursivity in the mathematical formulation of our al- 
gorithm given in Eq. (|3j). This means that we have to find 

a recursion formula for Pjy ' ■ • Pi |"0i,iv) 

\Si,...,S N ;rn N ) 

As shown in the previous example, only the fc'th sum- 
mand in Eq. ^ is unequal to zero, because the algorithm 
tells us to set up exactly k! = ^ + m N c~ polarizers 
(what means that exactly k' /3-coefficients are unequal 
zero and therefore the sum over k in Eq. |3]) reduces to 
k = k'). 

In the following we will divide Eq. ([3| into two sums. 
The first sum corresponds to Eq. ^ for N — 1 atoms 
times the projection of the iVth atom to the state |+) 
and the second sum corresponds to Eq. ^ for N — 1 
atoms times the projection of the iVth atom to the state 
) . Hereby, for the iVth atom, we can choose any of the 



' N 



± mjv) atoms which are projected into the state |±) 



iV 



± rn N ) 



thus we have to consider an additional factor ( 
for both sums. 

This decomposition can again be divided into N sums, 
since in principal every detector can project the iVth 
atom into the desired state. After this partitioning, we 
can exclude the factor fiiXiN {oi%XiN) what means that 
the iVth atom is projected into the state |+) (|— )) by 
registering the emitted photon at detector i. Thus, alto- 
gether we arrive at 



Pn ■ ■ ■ Pi \iPi,n} , — , v .. , , x 

fe=(f +m N ) (f +m N y.(f -m N )\ 

2J ' ' ' ^(f +m JV ) Q; ( T(f +m N + l) ■ ■ ■ a *(N)Xa(l)T(l) ' ' ' Xa(N)r(N) 



|+> T(1) (e| ® • • • ® |+) r(f +mjv) (e| ® |-) t( f+mjv+1 ) (e| ® • • • ® |-} r(JV) (e| J |^i,jv> 



G 



(f+mjy) 
(f + m^)!(f -mjv)! 



' ' ' Xa(JV)r(JV-l) (l+) T (i) ( e l ® ' ' ' ® l+) T (f + m „-l) (e| ® h) r (« +mjv ) (e| 
®---«|->t(jv-x) ( e l®l+M e l) lV>i,JV> 



+PnXnn X X /3 ( r(i)"-^(« +mjv _i)a (T (« +mjv )---a ( T(iV-i)X<T(i)r(i)--- 

• • • Xa(N-iMN-i) (l+)r(i) (e| <8> • • • ® |+) T ( f +mjv _i) (e| (8 h) T ( f +mjv ) (e| (10) 

®'"®|-)r(W-l) ( e l® l+>Ar(e|) IV'i.Jv) J 



' AT 
v 2 



(| + mjv )!(f - mjv )! 

aiXiw X X ^(2)---P^N +mN+l )a a ^N +mN+2 y--a <j{N) Xa(2)r(i)--- 

• • • X*{N)r(N-i) (l+)r(i) (e| ® • • • ® |+) r(f +mjv ) (e| ® |-) r(f +mjv+1 ) (e| 
®---®|->r(JV-i) h>7v(e|) IV'i.Jv) 



(11) 



+«ATXJVA' J! X/ ^(l)"-^(f+m N ) a <r(f +m JV + l) ' ' ' M»-l)X,r(l)r(l) ' ' ' 

• Xa(iv-i)x(iv-i) (l+)r(i) (e| ® • • • ® |+> T (f +mjv ) (e| ® |-> r ( f +mjv+l) (e 
® • ' ' ® l-)r(jv-i) (e| ® I-) at (e|) \ A.n) j • 



Hereby, we denoted with Sjv\(i) i n Ec^- (10) and (TOT) 
the symmetric group of the set {1, ...,N} and wiui 
Sat_! the symmetric group of the set 



{1, ...,N — 1}. Before simplifying Eqs. (10) and (11) fur- 
ther, we want to prove the following identity: 



When applying our algorithm in the mathematical formulation of Sect. Ill (Eq. (13])) to the generation of an arbitrary 
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N — 1 qubit total angular momentum eigenstate in an (N — l)-qubit compound system , the following equation holds: 



Pa(2) ■ ■ ■ Pv{k)&a(k+V) ' ' ' <V(jV-l)«<7(iV)X<7(2)T(l) ' ' ' Xct(N)t(N-1) 

crSSiv\(i) 

= 2^ ' ' ' P<r{k) a a{k+X) ' ' ' "<r(iV-l)Q!<7(JV)X(r(l)T(l) ' ' ' X<7(iV)T(iV-l) 



(12) 



#7(1) • • • /3 t r(fc-l) Q! <T()fe) ' ' ' Q ; ( t(A'-1)X<j(1)t(1) • • • Xo(N-1)t(N-1)- 



crGS 



N\(N) 



The identity follows from the prescription of the al- 
gorithm and the corresponding degrees of freedom: ac- 
cording to the algorithm, for the generation of the total 
angular momentum eigenstates for N — 1 qubits, we can 
freely choose any N — 1 out of the N detectors for the 



required setup. We exploit this degree of freedom by 
choosing a different detector configuration in each line of 
Eq.(12). In this way, we can substitute the JVth detector 
with any detector i £ {1, N — 1}. 

After that substitution we can also write 



#7(2) • • • Pa(k) a cj(k+1) ■ ■ ■ a C r(JV-l) Q! cr(l)X(T(2)r(l) ' ' ' Xa{l)r{N-l) 
y #7(1)#7(3) ■ ■ ■ ^<7(fc)Q!ff(fc+l) • • • Q!ff(JV-l)Q!<r(2)X<r(2)T(l) ' ' ' Xer(2)-r(jV- 

<7eSjv-i 



y #7(1) • • • Pa{k-l)®o{k) • • • Q!o-(JV-1)X<7(1)t(1) • • • Xo-(JV-1)t(JV-1) 



what proves Eq. ( 12 1. and in Eq. ( 11 ) (with k = y + + 1) we can simplify 

Eqs. ((lO) aoarljllj to: 

By using Eq. Q in Eq. (fToh (with fc = f + m N ) 



\k=(f +mN ) (f +m2V -l)!(f -m N )\ 

l^PiXiN] X /^(l)'' A(f+mW) a ff (f+m„)•■•^(^ r -l)^(lMl)••• 

Vi=l / reEjv-i ct6Ejv-i 

• • • Xa(N-i)r(N-i) (l+) T (i) (e| ® • • • ® | +) r ( f +mjv _i) (e| ® |-) r( f+mjv) (e| 

®"-®|-)r(tf-i) ( e l® l+M e l) IV'i.Jv) 



(f+ TOjv )!(f -m w -l)! 
[X)^^) X X / 3 o-(i)"-/ 3 (T (« +mjv )a (T (« +miv+ i)---a ( T(w-i)X CT (i)r(i)--- 

\i=l / tGSn-i ctGSjv-i 

• ' ' Xa(N-l)r(N-l) (l+)r(l) ( e l ® ' ' " ® l + )r(f +m N ) ( e l ® I ~> r ( f +m N + l) ( e l 



't(JV-I) 



We can reduce this expression further by using Eq. ^ for N - 1: 

Av • • A |V>i,JV> , N = ^AXiJvj Av-i—AlV^JV-i) 



fc=(f +mjv-l) 



1+) 



A' 



(13) 



+ X a *XiJV Av-1 ' ' ' A |^i,JV-l) 



According to Sect. II, in order to generate the state |Si, Sn] Wjv) we need k = {M- + to at) er —polarizers, for the 
state | Si, Sn-i', m N ~ h) we require k = +mjv — l) er polarizers and for the state Si, Sjv-i; Wat + |) 



we need fc = (tt + to a?) er —polarizers. Consequently, using Eq. (13 1, we can write the requested recursion formula: 



AT 



Pn - Pi \tl>i, 



\Si,...,S N ;mN) 



tew • Av-i • • ■ A |^i,iv-i) 



+ X a *XiJV ' Av-1 ' ' ' A IV'i./V-l) 



|jSi,...,5w_i;mw+|) 



1+) 



(14) 



r 



Next, we can compute the sums appearing in the 
brackets of Eq. ( 14 ) by considering the rules given by 



the algorithm in Sect. II. Hereby, we have to examine 
two cases: 



1. S N > S N _!: 

According to the algorithm, we have to connect 
the Nth atom with all detectors, if and only if 
these detectors weren't used for a descent Si < 



remaining coefficients as well as the dedicated 
phases \jk take all the values 1. Since k indicates 
the number of the descents, i.e., k — y — Sat 

Y 
2 



(Sat-i + |j, we can conclude that: 



JV 



JV 



'JV 
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Sj_i. As we have to set up exactly 
er polarizers and — 777at) <t+— polarizers, 
there are only (^+toat) /?— coefficients and 
(^r — mjv) a— coefficients unequal to zero. These 



»=i 



JV 



N 



+ mjv — k ) = [ Sat_i + 777at + 



(15) 



N 



SaT-1 - 777AT + - 



(16) 
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Sn < Sn-i: 

According to the algorithm we have to connect the 
N— th atom with two detectors, which are equipped 
with polarising filters of different orientation, i.e., 
one oriented along cr and the other along er+. 
Hence, only two phase factors XiN in our sums are 
unequal to zero. The connection between the iVth 
atom and the detector I which has a <r~ — polarizer 
in front (/?/ = 1 A ai — 0) must induce an 
optical phase shift of tt (xin = — 1)- If we de- 
note the other detector which is equipped with a 
<t + — polarizer in front with m (/3 m =0 A a m 



I A XmN — 1), we can conclude that: 

N 

^PiXiN = PlXlN + PmXmN = -1 (17) 
i=l 
N 

^ OtiXiN = OLlXlN + CtmXmN = 1 (18) 



With these identities we can formulate the requested 
recursion formula in the following way: 



Recursion formula for the algorithm: 

• Sn > Sn-i: 



• Sn < Sn-i: 



P n ■■■ Pi \i>i,N 

Sn-i + rriN 



\Si,...,S N \m N ) 
1 



P N -1 ■■■Pi \lpi,N-l) 



+ ( Sn-i - m N + ~ ) • Pn-i ---Pi \i>i,N-i 



|<Si,...,S«-_i;mjv — §) 
|Si,...,Sjv_i;-mjv+i) 



1 + ) 



(19) 



Pn---Pi \<Pi, 



\Si,...,S N ;m N ) 
P N -l---Pl\Tpi,N-l) 

P N -1 ■ ■ ■ A \4>i,N-l) 



|<5i,...,5jv"— i;mjv— §) 
|Si,...,Siv-i;mjv + 1) 



1 + ) 



(20) 



VI. RECURSIVITY OF THE COUPLING OF 
ANGULAR MOMENTA IN QUANTUM 
MECHANICS 



In quantum mechanics, the total angular momentum 
eigenstates | J, M) are written in the decoupled basis 
liij 32] mi, VH2) in the following form [20 



In this section we recapitulate the well-known proce- 
dure in quantum mechanics to obtain the eigenstates of a 
total angular momentum (coupled basis) from the eigen- 
states of the angular momenta to be coupled (bare or 
decoupled basis) and derive the corresponding recursion 
formulas. This will allow us in the next section to com- 
pare the two recursion formulas, the one from our algo- 
rithm (Eqs. |l9|) and (20 1 from Sect. V above) and the 



n 31 



one obtained from the quantum mechanical addition of 
angular momenta (see Eqs. (26) and (27) below), and to 
prove that both lead to the same result. In this way, we 
show that the algorithm presented in [16] and the one 
used in quantum mechanics to construct the eigenstates 
of a total angular momentum using the addition of an- 
gular momenta of N spin-1/2 particles are identical. 



\J, M )= ^2 (ji, 32;m 1 , m 2 \ J,M)\j 1 ,j 2 ;m 1 ,m 2 ) 

mi=-j'i m 2 =-j 2 

( 21 ) 

where the elements of the transformation matrix 
(ii)j2; mi, m2\ J,M) are the familiar Clebsch-Gordan- 
Coefficients (CGC) 20J . If we consider the well known 
selection rules 



M = mi + 777-2 



and |ii - j 2 \< J < ji +32- 



(22) 



(23) 



for the CGC, we can simplify Eq. (21), since in our case 
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mi takes only the values ± \ : 



VII. PROOF BY INDUCTION 



.n± 2 ,M 



Ji, 



M- 
1 



1 1 



1 1 1 

ji, 2 ;M- 2 , 2 



:M • 



1 1 

2'~2 



• 1 fcfj. 1 1 



(24) 



with J = ji ± |. In the case ji = | with A; € No and 
j2 = 5, the CGC are given by the following analytic 
expressions [20] : 



■ 1 nr 11 

n,-;M — , - 
J '2 2 2 



h,- 2 ,M- 



1 1 

2'~2 



■ 1 M 11 



1 1 1 

^'2 ;M+ 2'-2 



h + - 2 ,M 



.n + - r M 



h-- r M 



h~- 2 ,M 



2ji 



/Ji-M- 



2ii 



'Jl 



1 



2ji 



/ ji + j 
2ji + 1 



(25) 



If we denote the states |;±|) with |±) and recall that 
in our notation we have j\ — S'jv-i and j'2 = 5, we can 
write Eqs. (24) also in the following form: 



Recursion formula for the total angular momentum 
eigenstates with ji = Sjsr-i and j2 = | 

• Sn > Sn-i: 



By using the recursion formula for our algorithm 
(Eqs. (19) and p0| ) and the recursion formula for the 
coupled basis according to the quantum mechanical pre- 



scription for the addition of angular momenta (Eqs. (26) 



and (27)), we can now demonstrate the equivalence of 



the two expressions and in this way prove the validity of 
our algorithm by mathematical induction. 

Proposition: 

The algorithm presented in [16) is able to generate any 
of the 2 N symmetric and nonsymmetric total angular mo- 
mentum eigenstates spanning the Hilbert space of an ./V 
qubit compound, that is: 



\Si, S N ; m N ) = A - P N - - Px \%I)%,n) 

with AeR. 
Proof: 



|Si,...,S N ;mjv) 

(28) 



We prove this proposition with a mathematical 
induction for N . 



Basis: 



N = 2 



For N = 2 we can directly compute the proposition. 
We verify the proposition by using the well known states 
of the spin-1 triplet and the state of the spin-0 sin- 



and the rules of the 
with N = 2. 



gulet of a two qubit system [ 
algorithm as formulated in Eq. ([3 

As we can sec in the following table, we find for every 
state of a two qubit system one necessary ieK: 



I Si, Sn] itln) — 

'Sjv-i + m N + \ 
25at_i + 1 



Si, S N -i;m N - I ) <g) |+) 



Sn-i — rn,N + \ 
2S N - 1 + 1 



Si, Sn-i'^n + 



(26) 



• S N < S N _!: 

I Si, Sat; tojv) 



Sjv_i — 


m N + 


1 

2 


2Sn- 


-i + l 




Sjv-i + 


m N + 


1 

2 


2Sn- 


-1 + 1 





Si,...,S N -i;m N - - ) (g) |+) 



S , i,...,S , Ar_i;m Ar + -)<S)\-). 

(27) 
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Si, Sa; wjv) 


AA 1^,2) 

1 |Si,S 2 ;mjv> 


,4 


|a» L > -V 


I++; 


O 1 1 1 \ 

* P — w 


1 

2 


la> 


1 \ 1 1 1 w 
7i(l+-> + l-+» 


/1 1 \ 1 1 1 ^ 

U+-) + 1 — n) 


1 

^2 




l~> 


2|— ) 


1 

2 


|i,0;0> 


7S (!+->-!-+» 


(!+-)-!-+» 


1 

^2 



TABLE I: Verification of tie proposition with N — 2. 



r 



Induction hypothesis: |Si, <Sjv; wjv) by using Eqs. (26) and (27) of Sect. VI 

Let us assume that the following identity holds for one for the quantum mechanical procedure to generate the 



N £ N. 

\Si, Sjv-i; m N-i) = 

A - P N _i---Pi \ipi, N -i) 



|Si,...,Sjv-i;mjv_i) 



Induction step: 



N- 1 — >• TV 



Since the induction hypothesis holds for both states of 
the N — 1 qubit system: 



1 



Si, ■ Sjy^i;m N — - 



Ai ■ P N -i---Pi \lpi,N-l) 



\Si,...,Sn -i;mN— ^) 

(29) 



and 



|Si,...,Sjv-i;mjv + 5) 

(30) 



total angular momenta eigenstates. 

In the case of Sn > Sn-i, the state is given by (see 
Eq. d26b) 



Sn', tun) — 

Sn-i + m N + \ 
V 25jv-i + 1 



Si, ...,SN-i)mN - g ) ® 1+) 



Sjv-i — m w + 9 



Si,...,S N -i;rriN + - ) 



V 25at_i + 1 
whereas in the case of Sat < SjV-i the state can be ex- 



pressed by (see Eq. (27)) 
|5i, Sat; mjv) = 



5at_i — 


TOJV + 


1 

2 


2Sn- 


-1 + 1 




Sn-i + 


TOat + 


1 

2 


2i5jv- 


-1 + 1 





Si,..., S N -i;m N - - ) ® |+) 



Si,...,S N _i;m N + - ) (g) |-). 



If we now use the induction hypothesis (Eq. (29) 
with A\,Ai £ K, we can compute recursively the state and (30)), we obtain for Sn > Sjv-i 



' SjV_i + TOat + 1 

'SV-i - + I 



25jv_i + 1 



Ai • P 



jv-i • • • 



Pi 1^. 



i,N-l, 



25 



iV-1 



|Si,...,Sjv— i;mjv— 
A 2 ■ Pjv-i ■•■A |^,AT-l) . 

|Si,...,Siv-i;mjv + 5) 



1+) 



(31) 



I-) 
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and for Sn < Sn-x 

| Si, Sn', tun) 



' Sn-i - m N + | 
2SW-1 + 1 



Ai • P N -X ■■ - Pi \lpi,N-l 



2SW-1 + 1 



N- 



Next, we show that the right hand sides of Eqs. ( 
and (32 1 are, up to a real constant, nothing else t 



31) 



the right hand sides of Eqs. (19) and (20) corresponding 



to the recursion formula of our algorithm (see Sect. V). 
This identity proves the proposition Eq. ([28]). To this 
end, we insert Eq. (31) and ( [32] ) and the two recursion 
formula of the algorithm for Sn > Sn-x (Eq. (19)) and 



for S N < Sn-i (Eq. (|20j)) in the proposition Eq. (|28|. 
By comparing the coefficients of the appearing states, 
the proposition will be reduced to the following system 
of equations: 



For Sn > Sn-i the proposition reduces to 



Sn- 



771 JV 



2S N -i + 1 



A 1 = A[ S N - 



m N 



(33) 



" A ,_, " A+I .l; = .l|.S\ , n,, O (34) 



V 2S A r_ 1 + l 

and for Sn < Sn-x to 



' S N -i - m N + 5 
25iv_i + 1 



' Sn-i + m N + \ 
25jv-i + 1 



• A\ = -A 



A, = A. 



(35) 



(36) 



with A,AeR 
Special case: 

Sn-i + m N + \ = V Sn-i — tun + \ = 

In this special case, one of the two CGC appearing 
in the quantum mechanical recursion formula for state 
|jSi,...,Sjv;tojv) (Eqs. (26) and (27)) is zero. In that 
case the proposition Eq. (28) is obviously right since 



every equation in the form a ■ x = y (a ^ 0) has a unique 
solution in the field of the real numbers. In fact, we can 
easily indicate the requested A, A e R: 



|Si,...,5'iv-i; rn jv — h) 

1 ' ' ' A \lpi,N-l) . . , 

|Si,...,Sjv-i;w»jv+f ) 



(32) 



r 



If Sn-i + wjv + h = 0, A is given in the case of Sn > 
Sn-i by: 



A M '\j {2S N -x + 1) (Sn-! - m N + |) 
whereas for the case Sat < Sn-i, A is given by 



2S N -i + 1 



If Sn-i — mN + \ = 0, A is given in the case of Sat > 
Sn-i by 



A = Ax 



(2S N -i + 1) (SV-i + m N + |) ' 
whereas for Sn < Sn-i A is given by 



A = A 2 - 



' S N -i + m N + I 



2S N -i + l 
General case: 

Sn-i + m N + \ 7^ A Sjv-i — toat + \ ^ 

For all other cases, that is when all needed CGC are 
unequal to zero, we prove the proposition in the case of 
S N > Sn-i, if (see Eqs. Q and pi): 



1 



(2SjV-l + 1) (SaT-1 + 777AT + i) 



a?i := Ax 
is equal to 

^ '' = M ' ^ (2S N - 1 + l)(S N -i-m N + ^' 
In the case of Sn < Sn-i, we prove the proposition, if 



(see Eqs. (35) and (36)) 



&i = A\ 



'Sn-x - m N + \ 
2S N -x + 1 
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is equal to 



&2 — A 2 ■ 
i.e. in summary in all cases if 



S N -i + m N + 5 
2SW-i + 1 



holds. 



Ax _ S N -x +m N + | 
A2 V Sn-i - m N + \ 



In the following we retrieve this relation for Ax and 
A2 by employing a case differentiaton Sn-i > Sjv-2 and 

Sn-1 < 



First case: Sn—x > SW— 2 



(37) 

We start to use the recursion formula (Eq. ( |26| ) 
and (19l) for the states \Sx, Sjv-i ; Wjv — §) and 

Pjv-i • • - A |V>i,JV-i) . 

|5i,... 1 5 , Ar-i;miv — 2/ 



We thus obtain 



5i, Sn-x; tun — 



as well as 

P N -i ■••Pi iV'i.w-i) 



, 5 , AT_ 2 + mjv 
25 w _ 2 + 1 



|5i,...,S*Ar_ 2 ;mjv - 1) (8 |+) 



'/SjV-2 - + 1 

2S N - 2 + 1 



|S , l,...,5Ar_ 2 ;TOw) ® |-) 



|5'i,...,5'iv-i;"iiv — 



: (SjV-2 + m N ) Pjv-2 ' ' ' Pi \i>i,N-2) 

|Si,...,Sjv_2;mjv— 1) 
+ (Sjv-2 - niAr + 1) P w _ 2 • • • P x l^jv-2) 

|Si,...,Sjv-2;-mjv> 



I-) 



By using the induction hypothesis (Eq. ([29]) ) we can conclude from Eqs. (38) and (39) 



Sn-2 + m N . 1 1 \ 1 / SV-2 - tojv + 1 1 „ c , 1 . 

|6i, ...,S N - 2 ;m N - 1) ® |+) + W — — — — |6i, ...,b N -2;m N ) ® |-) 



y 25jv-2 + 1 

^1 ' (SV-2 + "Ijv) Av-a '"A IV'i.JV-a) 



25 w _ 2 + 1 

1+) 



|Si,...,S -N _ 2 ;mj V — 1) 
+ ^1 ' (SjV-2 - m N + 1) Av-2 • • • A \lpi,N-2> 



\S 1 ,...,S N - 2 ;m N ) 



r 



(38) 



(39) 



This equation implies that with 

the states \Sx, • SW-2; it^n) and 

Av-2 ■ ■ • A |"0»,iv-2) are linearly de- Pi 

|Si,...,Sjv-2;mw> 

pendent, i.e. it exists one pi £ K, so that the following 
equation holds: 



\Sx,—,S N - 2 ;m N ) =px-Pj 



N-2 ■ 



■ A 1^ 



i,JV— 2/ 



= Ai • V(2^- 2 + l)(^2-m N + l). (41) 



|Si,...,Siv-2;mjtf) 



In the same way we use the recursion for- 
mula for the states |iSi, Sn-x', tjjv + §) and 

A\r-i"-A|V»i,JV-i} , ( scc E qs- H 

Si,...,Sjv_i;mjv + 5 ) 



(40) and (19)). Hereby, we obtain 



Si, Sjsr-x;rriN + 



1\ Sn-2 + m N + 1 



2Sn-2 + 



|5i,...,S , Ar_ 2 ;"l A r) (gl |+) 



' jjV-2 - ^JV 
25^-2 + 1 



\Sx,...,S N - 2 ;m N + l)® |-) 
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as well as 



P 



V-i ■■■Pi \i>i,N-i) 



|Si,...,S'jv_i;mjv+|) 



[S N - 2 + m N + 1) P/V— 2 • • • A |^,AT- 2 ) ® |+) 

|Si,...,S^_ 2 ;rnjtf) 

+ (SjV-2 - Wat) Av-2 ■ ' ' A l^i.JV— 2) 
By using the induction hypothesis (Eq. (30)) we conclude that: 



|Si,...,Sjv-2;miv + l) 



/ Sn-2 + m N + 1 



2Sn-2 + 
^4 2 



- — \S 1 ,...,S N - 2 ;m N ) (8 |+) + 4 / 2g +1 
r_ 2 + mjv + 1) Av-2 ■ • • A 1^,^-2) 



|5i, S N - 2 ;m N + 1) 



(SW-2 + ™iV + 1) Av-2 ■ • • A \^i.N-2 

+ A 2 ■ (S N - 2 - m N ) P 



® 1+) 

|Si,...,Sjv-2;mjv> 

'v-2 • • • A IV'i.JV-a) 



|Si ,...,Sjv-2;mjv + l) 



r 



This equation again implies that 

the states \S\, Sn-2', tun) and 

P/v-i ' ' ' Pi |^i,JV- 2 ) are linearly dc- 

|Si,...,Sjv-2;mjv> 
pendent, i.e. it exists one p 2 € K, for which holds: 



and thus 



\Si,...,S N - 2 ;m N ) = p 2 - Pv-2 • • • Pi |V^,w-2> 
with 



|Si,...,Sjv_2;mjv) 
(42) 



1^2 



= A 2 ■ ^/{2S N - 2 + l)(S N „ 2 +m N + l). (43) 



Note that we can write the state \Si, Sn-2', tun) 
always as a linear combination of the decoupled basis 
k) (i€{l,...,2^- 2 }): 



\Si, S N ~ 2 ;m N ) = d\ 



with q € K. Beyond that, we can also write the state 



Ptv-2 ■ ■ • Pi lV^iV-2 
tion of the bare basis 



|Si,...,SjV-2,'TOjv) 



as a linear combina- 



Pv-2 ••Pi IVV-2) 



|Si,...,Sjv-2;mjv) 



with c?i e E. As we can write one state in the Hilbert 
space for a given basis in just one unique way, we can 
conclude that due to Eqs. (40) and (42) we have: 



Pi = -j- =P2- 



Hence, we obtain by using Eq. (41) and (43) 



Ayy/(2S N - 2 + 1) {S N _ 2 - m N + 1) = 

A 2 ■ ^(2S N -2 + 1) (Sn-2 + m N + l) 



Ai _ Sjy- 2 +rriN + l 
A 2 y Sn-2 - m N + 1 

As we assumed in this case Sjv-i > Sn-2, i.e. 

1 



S 



N-2 — Sn-1 



we have proven the equation 



A 2 



'Sn-i + m N + | 

Sn-1 — m N + 2 



(44) 



(37) 



for this case. 

Second case: Sn-i < Sn-2 

We can handle this case in the same way we 
treated the case Sn-i > Sn-2- By using the re- 
cursion formula for the states \Si, S/v-i; mN — 5) 



and Pn-i ■■■Pi \ipi,N-i) 



Si ,. . . ,S n — 



for Sn-i < 



Sn-2 ( sce Eqs. (|27j) and (20)), we obtain: 
1' 



Si, Sn-i;itin — 



'Sn-2 - m N + 1 
2S N -2 + 1 



\S 1 ,...,S N - 2 ;m N -l)(E)\+) 



S N -2 + m N , „ „ , , , 

+ 1 / -ttt; — r 1*1) ^at) |-) 



2Sn- 



as well as 
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Ptf-l • • ■ P WMr-l} , , , = - Pn-2 ■ ■ ■ Pi \i>i,N-2) ®|+) 

|5i,...,Sjv-i;mjv-5) |Si,...,Sjv-2;-tojv— 1) 



i.N-2, 



\Si,...,SN~2;m N ) 



I-) 



By means of the induction hypothesis Eq. (291 we can thus conclude that: 



'Sn-2 ~ m N + 1 I o i\o>li \ i SN-2 + m N i i \ 

|5i,...,5jv_ 2 ;mjv - 1) <8> |+) + 4/^7, — |-Si, -, S N -2]m N ) ® |-) 



25at_ 2 + 1 

- Ai ■ Pv_2 • • • Pi |*0i,JV- 2 ; 



|Si,...,Siv-2;mjv-l> 



25jv- 2 + 1 
8|+)+Ai- Pjv-2---Pi|Vi,iv-2) 



|Si,...,SAr_2;mjv) 



r 



This equation implies that the states l^i, S , at_ 2 ; m^) Pjv-i • • • A |^,,jv-i) in the case 

and Pn - 2 ■ ■ ' A '^^L,...^-,^)^ Unearly de " of S N -x < S N -2 (seeEqs'.' © and + |o|), we obtain: 
pendent, i.e. it exists one p\ € K, for which the following 
equation holds: 

1' 

\Si, ...,S N - 2 ;m N ) = pi-Pn-2 ■■■Pi \ipi,N-2) 

|Si,...,Sjv-2;mjv> 

(45) 



Si, S'at-i; totv + 



with 



' Sn-2 - m N 
2S N - 2 + 1 



\Si, S N -2;m N ) (g> |+) 



Pi = Ax 



' 2Sjy_ 2 + 1 
Sn-2 + rn N 



(46) 



'Sn-2 + m N + 1 
2S N ~2 + 1 



|S , i,...,5 Ar _ 2 ;m A r + 1) ® |-) 



Using again the recursion formula for 



the 



states 



and as well as 



Pv-i — PlV'i.Ar-i} . , = - Pn-2 ■ ■ ■ Pi \i>i,N-2) ®|+) 

|jSi,...,5^_i;mw+f ) |Si ,...,S«-2;"ijv) 



+ P W -2 ■••Pi IV'i.JV-a) 



|Si,...,Siv-2;niN + l> 



I-) 



By means of the induction hypothesis Eq. ( 30 ) we know that 



/S X 1 '" V S N . 2 ;m N ) |+) + + ^ + 1 \ Su Sjv _ a;row + l>® |-) 



2SjV-2 + 1 



- A 2 - Pv^-.-Pxl^ 



i,N-2i 



2S N -2 + 1 



|Si,...,Sjv-2;mjv> 



+ A 2 • Pjy-2 • • • Pi IV'i.AT-a) 



|Si,...,S , N -2;mjv + l) 



holds. 



Hereby, it is again implied that 
the states \S\, Sn-2] win) an d 



pendent, i.e. it exists one p 2 G R, for which holds: 



Ptv-2 • • • Pi \^i,N-i 



\Si,...,Snr-2;mnr) 



are linearly de- 



\Si, ...,Sn-2;itin) = P2-Pn-2 ■••Pi (^,^-2) 



|Sl,...,SjV-2;TOiv) 

(47) 
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with 



P2 = M 



2S N -2 + 1 



Sn- 



m N 



(48) 



As we can write one state in the Hilbert space for a given 
basis just in one unique way, we can conclude in the same 
way as in the first case that: 

Pi =P2- 



As we consider the case Sn-i < Sn-2, i.e. 

Sn-2 = Sjv-i + 2' 
we have again proven the equation 



Ai / Sn-i + rn N + \ 



Sn-i " m N + 



(37) 



Therewith, we obtain by means of Eqs. (46) and (48): 



A x _ S N -2 +m N 
A 2 y S N -2 - mjv 



The proven constraint for Ai and A 2 (Eq. (37)) is suf- 



ficient for the conclusion that a\ = a 2 as well as d\ = d 2 
holds. This means that altogether we can conclude that: 



ai = Ai 



{2S N -i + 1) (S N -i + m N + |) 



A 2 



! S N -i + m N + \ 



S N -i - m N + | Y (2S N -i + 1) (Sn-i + mN + |) 



=A 1 by using Eq. ]37[ 



= A 2 - 



(25jv_i + 1) {Sn-i -m N + \) 



= «2 



as well as that: 



di = A t 



A 9 



A 



2 ' 



'JV-1 



m N 



25jv_i + 1 



Sn-i + rnjsr + \ Sn-i 



S 



m N 



N-l 



m N 



2S N -i + 1 



by using Eq. \37\ 

Sn-i + m N + 5 
2Sjv_i + 1 



a 2 . 



Hence, the proposition is proven with A = oi\ = a 2 in 
the case of Sn > Sn-i & n d with A = d\ — d 2 in the 
case of Sn < Sn-i, respectively. 



tors such that the individual source of a registered pho- 
ton remains unknown. For the proof, we formulated the 
algorithm in form of a recursion formula what allowed 
us to compare it to the well-known quantum mechanical 
procedure for the generation of the eigenstates of a total 
angular momentum obtained by coupling a spin Sat_i 
with the spin of a further qubit. By comparing the two 
recursion formulas and demonstrating their equivalence, 
we verified the general validity of our scheme. A partic- 
ular feature of the algorithm is that it allows to encode 
the total spin eigenstates among independent and poten- 
tially far distant qubits. This means that the method is 
able to implement the required CGCs, needed to form 
those states, even though the qubits never directly inter- 
act with each other. In this way it simulates the coupling 
of the angular momenta of N non-interacting qubits. 



VIII. CONCLUSION AND OUTLOOK 

In this paper we have proven the algorithm introduced 
in [TB] enabling the generation of any of the 2 N symmetric 
and nonsymmetric total angular momentum eigenstates 
spanning the Hilbert space of an N qubit compound. The 
algorithm is implemented in a system of N single photon 
emitters with a A-configuration, scattering photons from 
an excited state |e) and decaying into ground states |+) 
and |— ), which form the individual qubit. Hereby it is 
assumed that the N photons are recorded by N detec- 



Note that in [5T] we presented a scheme for the gen- 
eration of all symmetric and non-symmetric total spin 
eigenstates encoded in the polarization degrees of free- 
dom of N photons. In the paper it was shown, that the 
methods for the generation of the polarization entangled 
photon states and for the generation of the entangled 
matter qubit states are mathematically equivalent. Thus, 
the proof of the algorithm given in this paper proves as 
well the algorithm derived in [21] for the generation of 
all symmetric and non-symmetric total spin eigenstates 
encoded in the polarization degrees of freedom of N pho- 
tons. 
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